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Abstract 

Two searches for supersymmetric particles in final states containing a same-flavour oppo¬ 
site-sign lepton pair, jets and large missing transverse momentum are presented. The proton- 
proton collision data used in these searches were collected at a centre-of-mass energy sfs - 8 
TeV by the ATLAS detector at the Large Hadron Collider and corresponds to an integrated 
luminosity of 20.3 fb -1 . Two leptonic production mechanisms are considered: decays of 
squarks and gluinos with Z bosons in the final state, resulting in a peak in the dilepton invari¬ 
ant mass distribution around the Z-boson mass; and decays of neutralinos (e.g. x\ —■ > f 

resulting in a kinematic endpoint in the dilepton invariant mass distribution. For the former, 
an excess of events above the expected Standard Model background is observed, with a signi¬ 
ficance of 3 standard deviations. In the latter case, the data are well-described by the expected 
Standard Model background. The results from each channel are interpreted in the context of 
several supersymmetric models involving the production of squarks and gluinos. 
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1. Introduction 


Supersymmetry (SUSY) [1-9] is an extension to the Standard Model (SM) that introduces supersymmetric 
particles (sparticles), which differ by half a unit of spin from their SM partners. The squarks (q) and 
sleptons ( l’) are the scalar partners of the quarks and leptons, and the gluinos (g) are the fermionic partners 
of the gluons. The charginos (xf with i = 1,2) and neutralinos (x° with i = 1,2,3,4) are the mass 
eigenstates (ordered from the lightest to the heaviest) formed from the linear superpositions of the SUSY 
partners of the Higgs and electroweak gauge bosons. SUSY models in which the gluino, higgsino and 
top squark masses are not much higher than the TeV scale can provide a solution to the SM hierarchy 
problem [10-15]. 

If strongly interacting sparticles have masses not higher than the TeV scale, they should be produced with 
observable rates at the Large Hadron Collider (LHC). In the minimal supersymmetric extension of the 
SM, such particles decay into jets, possibly leptons, and the lightest sparticle (LSP). If the LSP is stable 
due to R-parity conservation [15-19] and only weakly interacting, it escapes detection, leading to missing 
transverse momentum (p™ ss and its magnitude C | mss ) in the final state. In this scenario, the LSP is a 
dark-matter candidate [20,21], 

Leptons may be produced in the cascade decays of squarks and gluinos via several mechanisms. Here two 
scenarios that always produce leptons (electrons or muons) in same-flavour opposite-sign (SFOS) pairs 
are considered: the leptonic decay of a Z boson, Z —> £ + £~, and the decay x^ —> f + £~x {) v which includes 
contributions from x° —> —> ( + C~x { \ and^ —> Z*x\ —> ( + C~x {) v In models with generalised gauge- 

mediated (GGM) supersymmetry breaking with a gravitino LSP (G), Z bosons may be produced via the 
decay x° { —■► ZG. Z bosons may also result from the decay x® —> Zx° { , although the GGM interpretation 
with the decay x® —> ZG is the focus of the Z boson final-state channels studied here. The x® particle may 
itself be produced in the decays of the squarks or gluinos, e.g. q —» qx® and g —> qqx®- 

These two SFOS lepton production modes are distinguished by their - distributions of dilepton invariant 
mass (mu). The decay Z —> £ + t~ leads to a peak in the mu distribution around the Z boson mass, 
while the decay x® £ + £~X\ leads to a rising distribution in mu that terminates at a kinematic endpoint 
(“edge”) [22], because events with larger mu values would violate energy conservation in the decay of the 
j >2 particle. In this paper, two searches are performed that separately target these two signatures. A search 
for events with a SFOS lepton pair consistent with originating from the decay of a Z boson (on-Z search) 
targets SUSY models with Z boson production. A search for events with a SFOS lepton pair inconsistent 
with Z boson decay (off-Z search) targets the decay X® 

Previous searches for physics beyond the Standard Model (BSM) in the Z + jets + £™ ss final state have 
been performed by the CMS Collaboration [23,24]. Searches for a dilepton mass edge have also been 
performed by the CMS Collaboration [24,25]. In the CMS analysis performed with yfs = 8 TeV data 
reported in Ref. [24], an excess of events above the SM background with a significance of 2.6 standard 
deviations was observed. 

In this paper, the analysis is performed on the full 2012 ATLAS [26] dataset at a centre-of-mass energy of 
8 TeV, corresponding to an integrated luminosity of 20.3 fb 1 . 
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2. The ATLAS detector 


ATLAS is a multi-purpose detector consisting of a tracking system, electromagnetic and hadronic calor¬ 
imeters and a muon system. The tracking system comprises an inner detector (ID) immersed in a 2 T 
axial field supplied by the central solenoid magnet surrounding it. This sub-detector provides position and 
momentum measurements of charged particles over the pseudorapidity 1 range |//| < 2.5. The electromag¬ 
netic calorimetry is provided by liquid argon (LAr) sampling calorimeters using lead absorbers, covering 
the central region (|//| < 3.2). Hadronic calorimeters in the barrel region (|//| < 1.7) use scintillator tiles 
with steel absorbers, while the pseudorapidity range 1.5 < |//| < 4.9 is covered using LAr technology 
with copper or tungsten absorbers. The muon spectrometer (MS) has coverage up to |t/| < 2.7 and is built 
around the three superconducting toroid magnet systems. The MS uses various technologies to provide 
muon tracking and identification as well as dedicated muon triggering for the range \r/\ < 2.4. 

The trigger system [27] comprises three levels. The first of these (LI) is a hardware-based trigger that uses 
only a subset of calorimeter and muon system information. Following this, both the second level (L2) and 
event filter (EF) triggers, constituting the software-based high-level trigger, include fully reconstructed 
event information to identify objects. At L2, only the regions of interest in rj-(f> identified at LI are 
scrutinised, whereas complete event information from all detector sub-systems is available at the EF. 


3. Data and Monte Carlo samples 

The data used in this analysis were collected by ATLAS during 2012. Following requirements based on 
beam and detector conditions and data quality, the complete dataset corresponds to an integrated lumin¬ 
osity of 20.3 ftr 1 , with an associated uncertainty of 2.8 %. The uncertainty is derived following the same 
methodology as that detailed in Ref. [28], 

Dedicated high-transverse-momentum (pj) single-lepton triggers are used in conjunction with the lower- 
pj dilepton triggers to increase the trigger efficiency at high lepton p-\ . The required leading-lepton p-\ 
threshold is 25 GeV, whereas the sub-leading lepton threshold can be as low as 10 GeV, depending on 
the lepton p-\ threshold of the trigger responsible for accepting the event. To provide an estimate of 
the efficiency for the lepton selections used in these analyses, trigger efficiencies are calculated using tt 
Monte Carlo (MC) simulated event samples for leptons with p-\ > 14 GeV. For events where both leptons 
are in the barrel (endcaps), the total efficiency of the trigger configuration for a two-lepton selection is 
approximately 96 %, 88 % and 80 % (91 %, 92 % and 82 %) for ee, ep and pp events, respectively. 
Although the searches in this paper probe only same-flavour final states for evidence of SUSY, the ep 
channel is used to select control samples in data for background estimation purposes. 

Simulated event samples are used to validate the analysis techniques and aid in the estimation of SM back¬ 
grounds, as well as to provide predictions for BSM signal processes. The SM background samples [29-40] 

1 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector 
and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the t/-axis points 
upward. Cylindrical coordinates (r, <p) are used in the transverse plane, $ being the azimuthal angle around the beam pipe. The 

pseudorapidity is defined in terms of the polar angle 8 as r] = - In tan(0/2). The opening angle A R in 77 —0 space is defined as 
A R= V( A 7 ) 2 + ( A 0) 2 - 
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used are listed in Table 1, as are the parton distribution function (PDF) set, underlying-event tune and 
cross-section calculation order in a s used to normalise the event yields for these samples. Samples gener¬ 
ated with MadGraph5 1.3.28 [41] are interfaced with Pythia 6.426 [42] to simulate the parton shower. All 
samples generated using Powheg [43-45] use Pythia to simulate the parton shower, with the exception of 
the diboson samples, which use Pythia8 [46]. Sherpa [47] simulated samples use Sherpa’s own internal 
parton shower and fragmentation methods, as well as the Sherpa default underlying-event tune [47] . The 
standard ATLAS underlying-event tune, AUET2 [48], is used for all other samples with the exception of 
the Powheg-i-Pythia samples, which use the Perugia201 1C [49] tune. 

The signal models considered include simplified models and a generalised gauge-mediated supersymmetry¬ 
breaking model. In the simplified models, squarks and gluinos are directly pair-produced, and these sub¬ 
sequently decay to the LSP via two sets of intermediate particles. The squarks and gluinos decay with 
equal probability to the next-to-lightest neutralino or the lightest chargino, where the neutralino and char- 
gino arc mass-degenerate and have masses taken to be the average of the squark or gluino mass and 
the LSP mass. The intermediate chargino or neutralino then decays via sleptons (or sneutrinos) to two 
leptons of the same flavour and the lightest neutralino, which is assumed to be the LSP in these models. 
Here, the sleptons and sneutrinos are mass-degenerate and have masses taken to be the average of the 
chargino or neutralino and LSP masses. An example of one such process, pp —> gg —» (qqx^)(qqx^),- 
x\ —> t + t~x { \, xf —■> iG vx ( \ is illustrated on the left in Fig. 1, where £ = e,p,T with equal branch¬ 
ing fractions for each lepton flavour. The dilepton mass distribution for leptons produced from the x\ 
in these models is a rising distribution that terminates at a kinematic endpoint, whose value is given 
by nimax ~ mix 2 ) ~ m (Ti) = 1/2 (m(g/q) - mix®))- Therefore, signal models with small values of 
Am = m(g/q) - mix j) produce events with small dilepton masses; those with large Am produce events 
with large dilepton mass. 

For the model involving squark pair production, the left-handed partners of the u, d, c and s quarks have 
the same mass. The right-handed squarks and the partners of the b and t quarks are decoupled. For the 
gluino-pair model, an effective three-body decay for g —> qqx { \ is used, with equal branching fractions 
for q = u, d, c, s. Exclusion limits on these models are set based on the squark or gluino mass and 
the LSP mass, with all sparticles not directly involved in the considered decay chains effectively being 
decoupled. 

In the general gauge mediation models, the gravitino is the LSP and the next-to-lightest SUSY particle 
(NLSP) is a higgsino-like neutralino. The higgsino mass parameter, p, and the gluino mass are free 
parameters. The U(l) and SU(2) gaugino mass parameters, M\ and Mi, are fixed to be 1 TeV, and the 
masses of all other sparticles are set at ~ 1.5 TeV. In addition, p is set to be positive to make X] ZG 
the dominant NLSP decay. The branching fraction for^ —> ZG varies with tan/3, the ratio of the vacuum 
expectation value for the two Higgs doublets, and so two different values of tan /j are used. At tan (3 = 1.5, 
the branching fraction for x® —> ZG is large (about 97 %) [50], whereas setting tan /3 = 30 results in 
a considerable contribution (up to 40 %) from X® —» hG. In these models, li is the lightest CP-even 
SUSY Higgs boson, with m/, = 126 GeV and SM-like branching fractions. The dominant SUSY-particle 
production mode in these scenarios is the strong production of gluino pairs, which subsequently decay to 
the LSP via several intermediate particles. An example decay mode is shown in the diagram on the right 
in Fig. 1. The gravitino mass is set to be sufficiently small such that the NLSP decays are prompt. The 
decay length ctnlsp (where tnlsp is the lifetime of the NLSP) can vary depending on p, and is longest 
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at fi = 120 GeV, where it is 2 mm, decreasing to ctnlsp <0.1 mm for p > 150 GeV. The finite NLSP 
lifetime is taken into account in the MC signal acceptance and efficiency determination. 

All simplified models are produced using MadGraph5 1.3.33 with the CTEQ6L1 PDF set, interfaced with 
Pythia 6.426. The scale parameter for MLM matching [51] is set at a quarter of the mass of the lightest 
strongly produced sparticle in the matrix element. The SUSY mass spectra, gluino branching fractions and 
the gluino decay width for the GGM scenarios are calculated using Suspect 2.41 [52] and Sdecay 1.3 [53]. 
The GGM signal samples are generated using Pythia 6.423 with the MRST20O7 LO* [54] PDF set. The 
underlying event is modelled using the AUET2 tune for all signal samples. Signals are normalised to 
cross sections calculated at next-to-leading order (NLO) in a s , including the resummation of soft gluon 
emission at next-to-leading-logarithmic accuracy (NLO+NLL) [55-59]. 

A full ATLAS detector simulation [60] using GEANT4 [61] is performed for most of the SM background 
MC samples. The signal and remaining SM MC samples use a fast simulation [62], which employs a 
combination of a parameterisation of the response of the ATLAS electromagnetic and hadronic calorimet¬ 
ers and GEANT4. To simulate the effect of multiple pp interactions occurring during the same (in-time) 
or a nearby (out-of-time) bunch-crossing, called pile-up, minimum-bias interactions are generated and 
overlaid on top of the hard-scattering process. These are produced using Pythia8 with the A2 tune [63]. 
MC-to-data corrections are made to simulated samples to account for small differences in lepton identi¬ 
fication and reconstruction efficiencies, and the efficiency and misidentification rate associated with the 
algorithm used to distinguish jets containing 6 -hadrons. 
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Figure 1: Decay topologies for example signal processes. A simplified model involving gluino pair production, with 
the gluinos following two-step decays via sleptons to neutralino LSPs is shown on the left. The diagram on the right 
shows a GGM decay mode, where gluinos decay via neutralinos to gravitino LSPs. 


4. Physics object identification and selection 

Electron candidates are reconstructed using energy clusters in the electromagnetic calorimeter matched 
to ID tracks. Electrons used in this analysis are assigned either “baseline” or “signal” status. Baseline 
electrons are required to have transverse energy Ej > 10 GeV, satisfy the “medium” criteria described in 
Ref. [64] and reside within I 77 I < 2.47 and not in the range 1.37 < |//| < 1.52. Signal electrons are further 
required to be consistent with the primary vertex and isolated with respect to other objects in the event, 
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Physics process Generator Parton Cross section Tune PDF set 

Shower 


Z/y*(—» Ct) + jets 

Sherpa 1.4.1 

Sherpa 1.4.1 

NNLO [29,30] 

Sherpa default 

NLO CT10 [31] 

tt 

Powheg-Box r2129 

Pythia 6.426 

NNLO+NNLL [32,33] 

Perugia201 1C 

NLO CT10 

Single-top ( Wt ) 

Powheg-Box rl556 

Pythia 6.426 

Approx. NNLO [34,35] 

Perugia201 1C 

NLO CT10 

t + Z 

MadGraph5 1.3.28 

Pythia 6.426 

LO 

AUET2 

CTEQ6L1 [36] 

tt + W and tt + Z 

MadGraph5 1.3.28 

Pythia 6.426 

NLO [37,38] 

AUET2 

CTEQ6L1 

tt + ww 

MadGraph5 1.3.28 

Pythia 8.165 

LO 

AUET2 

CTEQ6L1 

ww, 

WZ and ZZ 

POWHEG-BOX rl508 

Pythia 8.163 

NLO [39,40] 

AUET2 

NLO CT10 


Table 1: Simulated background event samples used in this analysis with the corresponding generator, cross-section 
order in » s used to normalise the event yield, underlying-event tune and PDF set. 


with a p-p-dcpcndcnt isolation requirement. The primary vertex is defined as the reconstructed vertex 
with the highest £ p 2 , where the summation includes all particle tracks with pj > 400 MeV associated 
with a given reconstructed vertex. Signal electrons with Ej < 25 GeV must additionally satisfy the 
more stringent shower shape, track quality and matching requirements of the “tight” selection criteria in 
Ref. [64]. For electrons with Ej < 25 GeV (> 25 GeV), the sum of the transverse momenta of all charged- 
particle tracks with pj > 400 MeV associated with the primary vertex, excluding the electron track, within 
A R = 0.3 (0.2) surrounding the electron must be less than 16 % (10 %) of the electron pj. Electrons with 
Ej < 25 GeV must reside within a distance \zo sin 9\ < 0.4 mm of the primary vertex along the direction 
of the beamline 2 . The significance of the transverse-plane distance of closest approach of the electron to 
the primary vertex must be |<7o/cA/J < 5. For electrons with Ej > 25 GeV, |zol is required to be < 2 mm 
and k/ol < I mm. 

Baseline muons are reconstructed from either ID tracks matched to a muon segment in the muon spectro¬ 
meter or combined tracks formed both from the ID and muon spectrometer [65] . They are required to be 
of good quality, as described in Ref. [66], and to satisfy pj > 10 GeV and |//| < 2.4. Signal muons are 
further required to be isolated, with the scalar sum of the p \ of charged particle tracks associated with the 
primary vertex, excluding the muon track, within a cone of size A R < 0.3 surrounding the muon being 
less than 12 % of the muon pj for muons with pj < 25 GeV. For muons with pj > 25 GeV, the scalar sum 
of the pj of charged-particle tracks associated with the primary vertex, excluding the muon track, within 
A R < 0.2 surrounding the muon must be less than 1.8 GeV. Signal muons with p \ < 25 GeV must also 
have |zosinF| < 1 mm and \do/cr^ 0 \ < 3. For the leptons selected by this analysis, the do requirement is 
typically several times less restrictive than the Ir/o/cr^J requirement. 

Jets are reconstructed from topological clusters in the calorimeter using the anti-A) algorithm [67] with 
a distance parameter of 0.4. Each cluster is categorised as being electromagnetic or hadronic in origin 
according to its shape [68], so as to account for the differing calorimeter response for electrons/photons 
and hadrons. A cluster-level correction is then applied to electromagnetic and hadronic energy deposits 
using correction factors derived from both MC simulation and data. Jets are corrected for expected pile-up 
contributions [69] and further calibrated to account for the calorimeter response with respect to the true 


2 The distance of closest approach between a particle object and the primary vertex in the longitudinal (transverse) plane is 
denoted by zo (do)- 
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jet energy [70, 71]. A small residual correction is applied to the jets in data to account for differences 
between response in data and MC simulation. Baseline jets are selected with px > 20 GeV. Events in 
which these jets do not pass specific jet quality requirements are rejected so as to remove events affected 
by detector noise and non-collision backgrounds [72,73]. Signal jets are required to satisfy px > 35 GeV 
and |//| < 2.5. To reduce the impact of jets from pileup to a negligible level, jets with px < 50 GeV within 
\rj\ < 2.4 are further required to have a jet vertex fraction |JVF| > 0.25. Here the JVF is the px-weighted 
fraction of tracks matched to the jet that are associated with the primary vertex [74], with jets without any 
associated tracks being assigned JVF = -1. 

The MV1 neural network algorithm [75] identifies jets containing /?-hadrons using the impact parameters 
of associated tracks and any reconstructed secondary vertices. For this analysis, the working point corres¬ 
ponding to a 60 % efficiency for tagging Zr-jets in simulated tt events is used, resulting in a charm quark 
rejection factor of approximately 8 and a light quark/gluon jet rejection factor of about 600. To ensure that 
each physics object is counted only once, an overlap removal procedure is applied. If any two baseline 
electrons reside within A R = 0.05 of one another, the electron with lower Ej is discarded. Following this, 
any baseline jets within A R = 0.2 of a baseline electron are removed. After this, any baseline electron 
or muon residing within A R = 0.4 of a remaining baseline jet is discarded. Finally, to remove electrons 
originating from muon bremsstrahlung, any baseline electron within A R = 0.01 of any remaining baseline 
muon is removed from the event. 

The Z2| llss is defined as the magnitude of the vector sum of the transverse momenta of all photons, elec¬ 
trons, muons, baseline jets and an additional “soft term” [76]. The soft term includes clusters of energy 
in the calorimeter not associated with any calibrated object, which are corrected for material effects and 
the non-compensating nature of the calorimeter. Reconstructed photons used in the E™ ss calculation are 
required to satisfy the “tight” requirements of Ref. [77], 


5. Event selection 

Events selected for this analysis must have at least five tracks with px > 400 MeV associated with the 
primary vertex. Any event containing a baseline muon with |zo sin 8\ > 0.2 mm or |r/o| > 1 mm is rejected, 
to remove cosmic-ray events. To reject events with fake E™ ss , those containing poorly measured muon 
candidates, characterised by large uncertainties on the measured momentum, are also removed. If the 
invariant mass of the two leading leptons in the event is less than 15 GeV the event is vetoed to suppress 
low-mass particle decays and Drell-Yan production. 

Events are required to contain at least two signal leptons (electrons or muons). If more than two signal 
leptons are present, the two with the largest values of px are selected. These leptons must pass one of 
the leptonic triggers, with the two leading leptons being matched, within A R < 0.15, to the online trigger 
objects that triggered the event in the case of the dilepton triggers. For events selected by a single-lepton 
trigger, one of the two leading leptons must be matched to the online trigger object in the same way. 
The leading lepton in the event must have px > 25 GeV and the sub-leading lepton is required to have 
px > 10-14 GeV, depending on the px theshold of the trigger selecting the event. For the off-Z analysis, 
the sub-leading lepton px threshold is increased to 20 GeV. This is done to improve the accuracy of 
the method for estimating flavour-symmetric backgrounds, discussed in Sect. 6.2, in events with small 
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dilepton invariant mass. For the same reason, the m/c threshold is also raised to 20 GeV in this search 
channel. The two leading leptons must be oppositely charged, with the signal selection requiring that 
these be same-flavour (SF) lepton pairs. The different-flavour (DF) channel is also exploited to estimate 
certain backgrounds, such as that due to tt production. All events are further required to contain at least 
two signal jets, since this is the minimum expected jet multiplicity for the signal models considered in this 
analysis. 

Three types of region are used in the analysis. Control regions (CRs) are used to constrain the SM back¬ 
grounds. These backgrounds, estimated in the CRs, are first extrapolated to the validation regions (VRs) 
as a cross check and then to the signal regions (SRs), where an excess over the expected background is 
searched for. 

GGM scenarios are the target of the on-Z search, where the G from —> (Z/h) + G decays is expected to 
result in Z?™ ss . The Z boson mass window used for this search is 81 < mu < 101 GeV. To isolate GGM 
signals with high gluino mass and high jet activity the on-Z SR, SR-Z, is defined using requirements on 
Z?™ ss and Hj = p^ tJ + /x^ 1110 ' 1 ' 1 + /?^ pton ' 2 , where Hj includes all signal jets and the two leading leptons. 
Since b -jets are often, but not always, expected in GGM decay chains, no requirement is placed on b- 
tagged jet multiplicity. Dedicated CRs are defined in order to estimate the contribution of various SM 
backgrounds to the SR. These regions are constructed with selection criteria similar to those of the SR, 
differing either in mil or MET ranges, or in lepton flavour requirements. A comprehensive discussion of 
the various methods used to perform these estimates follows in Sect. 6. For the SR and CRs, detailed in 
Table 2, a further requirement on the azimuthal opening angle between each of the leading two jets and the 
Zi™ ss (A0(jct | 2 , /?.'| llss )) is introduced to reject events with jet mismeasurements contributing to large fake 
Ej lss . This requirement is applied in the SR and two CRs used in the on-Z search, all of which have high 
E'™ ss and Hj thresholds, at 225 GeV and 600 GeV, respectively. Additional VRs are defined at lower E™ ss 
and Hj to cross-check the SM background estimation methods. These are also sumarised in Table 2. The 
SR selection results in an acceptance times efficiency of 2M %, including leptonic Z branching fractions, 
for GGM signal models with p > 400 GeV. 

In the off-Z analysis, a search is performed in the Z boson sidebands. The Z boson mass window vetoed 
here is larger than that selected in the on-Z analysis (mu t [SO, 110] GeV) to maximise Z boson rejection. 
An asymmetric window is chosen to improve the suppression of boosted Z —> pp events with muons whose 
momenta are overestimated, leading to large £™ ss . In this search, four SRs are defined by requirements on 
jet multiplicity, ^-tagged jet multiplicity, and Z2.' [ niss . The SR requirements are optimised for the simplified 
models of pair production of squarks (requiring at least two jets) and gluinos (requiring at least four jets) 
discussed in Sect. 3. Two SRs with a b -veto provide the best sensitivity in the simplified models considered 
here, since the signal b -jet content is lower than that of the dominant (/"background. Orthogonal SRs with 
a requirement of at least one b-tagged jet target other signal models not explicitly considered here, such 
as those with bottom squarks that are lighter than the other squark flavours. For these four SRs, the 
requirement E ™ ss > 200 GeV is imposed. In addition, one signal region with requirements similar to 
those used in the CMS search [24] is defined (SR-loose). These SRs and their respective CRs, which have 
the same jet and /F | 11ISS requirements, but select different mu ranges or lepton flavour combinations, are 
defined in Table 3. 

The most sensitive off-Z SR for the squark-pair (gluino-pair) model is SR-2j-bveto (SR-4j-bveto). Because 
the value of the mu kinematic endpoint depends on unknown model parameters, the analysis is performed 


over multiple mu ranges for these two SRs. The dilepton mass windows considered for the SR-2j-bveto 
and SR-4j-bveto regions are presented in Sect. 9. For the combined ee + Hiu channels, the typical signal 
acceptance times efficiency values for the squark-pair (gluino-pair) model in the SR-2j-bveto (SR-4j- 
bveto) region are 0.1-10 % (0.1-8 %) over the full dilepton mass range. 

The on-Z and off-Z searches are optimised for different signal models and as such are defined with ortho¬ 
gonal SRs. Given the different signatures probed, there are cases where the CR of one search may overlap 
with the SR of the other. Data events that fall in the off-Z SRs can comprise up to 60 % of the top CR 
for the on-Z analysis (CRT, defined in Table 2). Data events in SR-Z comprise up to 36 % of the events 
in the CRs with 80 < mu < 110 GeV that are used to normalise the Z + jets background in the off-Z 
analysis, but the potential impact on the background prediction is small because the Z + jets contribution 
is a small fraction of the total background. For the following analysis, each search assumes only signal 
contamination from the specific signal model they are probing. 


On-Z 

Region 

jjmiss 

[GeV] 

H T 

[GeV] 

Jfjets 

nt(( 

[GeV] 

SF/DF 

Ef ss sig. 

[ V GeV] 

/ST 

A*(jet 12 ,F” iss ) 

Signal regions 

SR-Z 

> 225 

> 600 

> 2 

81 < m (e < 101 

SF 

- 

- 

> 0.4 

Control regions 

Seed region 

- 

> 600 

> 2 

81 < mu < 101 

SF 

< 0.9 

< 0.6 

- 

CR e/r 

> 225 

> 600 

> 2 

81 < m (e < 101 

DF 

- 

- 

> 0.4 

CRT 

> 225 

> 600 

> 2 

m tc i [81,101] 

SF 

- 

- 

> 0.4 

Validation regions 

VRZ 

< 150 

> 600 

> 2 

81 < mu < 101 

SF 

- 

- 

- 

VRT 

150-225 

> 500 

> 2 

m f ( i [81,101] 

SF 

- 

- 

> 0.4 

VRTZ 

150-225 

> 500 

> 2 

81 < mu < 101 

SF 

- 

- 

> 0.4 


Table 2: Overview of all signal, control and validation regions used in the on-Z search. More details are given in the 
text. The E™ ss significance and the soft-term fraction fsi needed in the seed regions for the jet smearing method 
are defined in Sect. 6.1. The flavour combination of the dilepton pair is denoted as either “SF” for same-flavour or 
“DF” for different flavour. 


6. Background estimation 

The dominant background processes in the signal regions, and those that are expected to be most difficult 
to model using MC simulation, are estimated using data-driven techniques. With SRs defined at large 
f?™ ss , any contribution from Z/y* + jets will be a consequence of artificially high C | n ' ss in the event due 
to, for example, jet mismeasurements. This background must be carefully estimated, particularly in the 
on-Z search, since the peaking Z/y* + jets background can mimic the signal. This background is expected 
to constitute, in general, less than 10 % of the total background in the off-Z SRs and have a negligible 
contribution to SR-Z. 
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Off-Z 

jjmiss 

^jets 

^b-jets 

mu 

SF/DF 

Region 

IGeVj 



[GeV] 


Signal regions 

SR-2j-bveto 

> 200 

> 2 

= 0 

mu i [80, 110] 

SF 

SR-2j-btag 

> 200 

> 2 

> i 

m u t [80,110] 

SF 

SR-4j-bveto 

> 200 

> 4 

= 0 

mu t [80,110] 

SF 

SR-4j-btag 

> 200 

> 4 

> 1 

m u t [80, 110] 

SF 

SR-loose 

> (150, 100) 

(2,2s 3) 

- 

m f( t [80, 110] 

SF 

Control regions 

CRZ-2j-bveto 

> 200 

> 2 

= 0 

80 < mu <110 

SF 

CRZ-2j-btag 

> 200 

> 2 

> 1 

80 < m u <110 

SF 

CRZ-4j-bveto 

> 200 

> 4 

= 0 

80 < m u <110 

SF 

CRZ-4j-btag 

> 200 

> 4 

> 1 

80 < mu <110 

SF 

CRZ-loose 

> (150,100) 

(2,2s 3) 

- 

80 < m u <110 

SF 

CRT-2j-bveto 

> 200 

> 2 

= 0 

m u t [80,110] 

DF 

CRT-2j-btag 

> 200 

> 2 

> 1 

% t [80,110] 

DF 

CRT-4j-bveto 

> 200 

> 4 

= 0 

tti(£ [80,110] 

DF 

CRT-4j-btag 

> 200 

> 4 

> 1 

m u t [80, 110] 

DF 

CRT-loose 

> (150, 100) 

(2, >3) 

- 

m u t [80, 110] 

DF 

Validation regions 

VR-offZ 

100-150 

= 2 

- 

mu t [80,110] 

SF 


Table 3: Overview of all signal, control and validation regions used in the ofF-Z analysis. For SR-loose, events with 
two jets (at least three jets) are required to satisfy £™ ss > 150 (100) GeV. Further details are the same as in Table 2. 
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In both the off-Z and on-Z signal regions, the dominant backgrounds come from so-called “flavour-sym¬ 
metric” processes, where the dileptonic branching fractions to ee , pp and ep have a 1:1:2 ratio such that 
the same-flavour contributions can be estimated using information from the different-flavour contribu¬ 
tion. This group of backgrounds is dominated by tt and also includes WW, single top (Wt) and Z —> rr 
production, and makes up ~ 60 % (~ 90 %) of the predicted background in the on-Z (off-Z) SRs. 

Diboson backgrounds with real Z boson production, while small in the off-Z regions, contribute up to 
25 % of the total background in the on-Z regions. These backgrounds are estimated using MC simulation, 
as are “rare top” backgrounds, including tt+ W{W)/Z (i.e. tt + W,tt + Z and tt + WW) and t + Z processes. 
All backgrounds that are estimated from MC simulation are subject to carefully assessed theoretical and 
experimental uncertainties. 

Other processes, including those that might be present due to mis-reconstructed jets entering as leptons, 
can contribute up to 10 % (6 %) in the on-Z (off-Z) SRs. The background estimation techniques followed 
in the on-Z and off-Z searches are similar, with a few well-motivated exceptions. 


6.1. Estimation of the Z/y*+ jets background 

6.1.1. Z/y* + jets background in the off-Z search 

In the off-Z signal regions, the background from Z/y* + jets is due to off-shell Z bosons and photons, or to 
on-shell Z bosons with lepton momenta that are mismeasured. The region with dilepton mass in the range 
80 < mu <110 GeV is not considered as a search region. To estimate the contribution from Z/y* + jets 
outside of this range, dilepton mass shape templates are derived from Z/y* + jets MC events. These shape 
templates are normalised to data in control regions with the same selection as the corresponding signal 
regions, but with the requirement on mu inverted to 80 < mu <110 GeV, to select a sample enriched in 
Z/y* + jets events. These CRs are defined in Table 3. 


6.1.2. Z/y* + jets background in the on-Z search 

The assessment of the peaking background due to Z/y* + jets in the on-Z signal regions requires careful 
consideration. The events that populate the signal regions result from mismeasurements of physics objects 
where, for example, one of the final-state jets has its energy underestimated, resulting in an overestimate 
of the total E ™ ss in the event. Due to the difficulties of modelling instrumental C | lllss in simulation, MC 
events are not relied upon alone for the estimation of the Z/y* + jets background. A data-driven technique 
is used as the nominal method for estimating this background. This technique confirms the expectation 
from MC simulation that the Z + jets background is negligible in the SR. 

The primary method used to model the Z/y* + jets background in SR-Z is the so-called “jet smearing” 
method, which is described in detail in Ref. [78]. This involves defining a region with Z/y* + jets events 
containing well-measured jets (at low E™ ss ), known as the “seed” region. The jets in these events are then 
smeared using functions that describe the detector’s jet pj response and (p resolution as a function of jet pj, 
creating a set of pseudo-data events. The jet-smearing method provides an estimate for the contribution 
from events containing both fake if™ 58 , from object mismeasurements, and real £™ ss , from neutrinos in 
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heavy-flavour quark decays, by using different response functions for light-flavour and ^-tagged jets. The 
response function is measured by comparing generator-level jet pj to reconstructed jet p \ in Pythia8 dijet 
MC events, generated using the CT10 NLO PDF set. This function is then tuned to data, based on a dijet 
balance analysis in which the pj asymmetry is used to constrain the width of the Gaussian core. The non- 
Gaussian tails of the response function are corrected based on > 3-jet events in data, selected such that the 
£™ ss in each event points either towards, or in the opposite direction to one of the jets. This ensures that 
one of the jets is clearly associated with the £™ 88 , and the jet response can then be described in terms of 
the E™ 88 and reconstructed jet p-\. This procedure results in a good estimate of the overall jet response. 

In order to calculate the if™ 88 distribution of the pseudo-data, the £™ 88 is recalculated using the new 
(smeared) jet pj and <p. The distribution of pseudo-data events is then normalised to data in the low- 
£™ ss re gj on (| () < £™ ss < 50 GeV) of a validation region, denoted VRZ, after the requirement of 
A(/)(jet[ 2 ,£™ ss ) > 0.4. This is defined in Table 2 and is designed to be representative of the signal re¬ 
gion but at lower E' | T1ISS , where the contamination for relevant GGM signal models is expected to be less 
than 1 %. 

The seed region must contain events with topologies similar to those expected in the signal region. To 
ensure that this is the case, the Hj and jet multiplicity requirements applied to the seed region remain 
the same as in the signal region, while the £™ 88 threshold of 225 GeV is removed, as shown in Table 2. 
Although the seed events should have little to no £™ ss , enforcing a direct upper limit on £™ 88 can intro¬ 
duce a bias in the jet p j distribution in the seed region compared with the signal region. To avoid this, a 
requirement on the £™ 88 significance, defined as: 


^■miss 

E? ss sig. = — T (1) 

4 el + £ Ef l 

is used in the seed region. Here X E^ 1 and X Ej 0ft are the summed Ej from the baseline jets and the 
low-energy calorimeter deposits not associated with final-state physics objects, respectively. Placing a 
requirement on this variable does not produce a shape difference between jet pj distributions in the seed 
and signal regions, while effectively selecting well-balanced Z/y* + jets events in the seed region. This 
requirement is also found to result in no event overlap between the seed region and SR-Z. 

In the seed region an additional requirement is placed on the soft-term fraction, /st, defined as the fraction 
of the total £™ 88 in an event originating from calorimeter energy deposits not associated with a calibrated 
lepton or jet (/st = X E' | msv Soft /E™ 88 ), to select events with small /st- This is useful because events with 
large values of fake £™ 88 tend to have low soft-term fractions (/st < 0.6). 

The requirements on the E™ 8S significance and /st are initially optimised by applying the jet smearing 
method to Z/y* + jets MC events and testing the agreement in the £™ 88 spectrum between direct and 
smeared MC events in the VRZ. This closure test is performed using the response function derived from 
MC simulation. 

The Z/y* + jets background predominantly comes from events where a single jet is grossly mismeasured, 
since the mismeasurement of additional jets is unlikely, and can lead to smearing that reduces the total 
£™ 1SS . The requirement on the opening angle in / between either of the leading two jets and the E™ ss , 
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AfWjet, 2 ,£'™ ss ) > 0.4, strongly suppresses this background. The estimate of the Z/y* + jets background 
is performed both with and without this requirement, in order to aid in the interpretation of the results 
in the SR, as described in Sect. 8. The optimisation of the Ej lss significance and /st requirements are 
performed separately with and without the requirement, although the optimal values are not found to 
differ significantly. 

The jet smearing method using the data-corrected jet response function is validated in VRZ, compar¬ 
ing smeared pseudo-data to data. The resulting £™ ISS distributions show agreement within uncertainties 
assessed based on varying the response function and the £! [ T1ISS significance requirement in the seed re¬ 
gion. The £' | mss distribution in VRZ, with the additional requirement A0(jet 1 2 , £™ ss ) > 0.4, is shown 
in Fig. 2. Flere the £' | mss range extends only up to 100 GeV, since tt events begin to dominate at higher 
£'| 1ISS values. The pseudo-data to data agreement in VRZ motivates the final determination of the /2' [ n ' ss 
significance requirement used for the seed region (Ej lss sig. < 0.9). Backgrounds containing real £'’ 1ISS , 
including tt and diboson production, are taken from MC simulation for this check. The chosen values are 
detailed in Table 2 with a summary of the kinematic requirements imposed on the seed and Z validation 
region. Extrapolating the jet smearing estimate to the signal regions yields the results detailed in Table 4. 
The data-driven estimate is compatible with the MC expectation that the Z + jets background contributes 
significantly less than one event in SR-Z. 


Signal region 

Jet-smearing 

Z+jets MC 

SR-Z ee 

SR-Z nn 

0.05 ± 0.04 

0 02 +0 - 03 

U - 0.02 

0.05 ± 0.03 

0.09 ± 0.05 


Table 4: Number of Z/y* +jets background events estimated in the on-Z signal region (SR-Z) using the jet smearing 
method. This is compared with the prediction from the Sherpa MC simulation. The quoted uncertainties include 
those due to statistical and systematic effects (see Sect. 7). 


6.2. Estimation of the flavour-symmetric backgrounds 

The dominant background in the signal regions is tt production, resulting in two leptons in the final state, 
with lesser contributors including the production of dibosons (WW), single top quarks ( Wt) and Z bosons 
that decay to r leptons. For these the so-called “flavour-symmetry” method can be used to estimate, in a 
data-driven way, the contribution from these processes in the same-flavour channels using their measured 
contribution to the different-flavour channels. 


6.2.1. Flavour-symmetric background in the on-Z search 

The flavour-symmetry method uses a control region, CRe/j in the case of the on-Z search, which is defined 
to be identical to the signal region, but in the different-flavour e/j channel. In CRe/u, the expected contam¬ 
ination due to GGM signal processes of interest is < 3 %. 
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Figure 2: Distribution of E™ ss in the electron (left) and muon (right) channel in VRZ of the on-Z analysis following 
the requirement of Ar/>(jet| 2 , £,’j nlss ) > 0.4. Here the Z/y* + jets background (solid blue) is modelled using pj- and 
0-smeared pseudo-data events. The hatched uncertainty band includes the statistical uncertainty on the simulated 
event samples and the systematic uncertainty on the jet-smearing estimate due to the jet response function and the 
seed selection. The backgrounds due to WZ , ZZ or rare top processes, as well as from lepton fakes, are included 
under “Other Backgrounds”. 
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The number of data events observed (N^ ata ) in this control region is corrected by subtracting the expected 
contribution from backgrounds that are not flavour symmetric. The background with the largest impact 
on this correction is that due to fake leptons, with the estimate provided by the matrix method, described 
in Sect. 6.3, being used in the subtraction. All other contributions, which include WZ, ZZ, iZ and it 
+W(W)/Z processes, are taken directly from MC simulation. This corrected number, V d /f a ’ corr , is related 
to the expected number in the same-flavour channels, , by the following relations: 


n: 


est 


n: 




_j_ »rdata,corr 
2 

_j_ Ardata.corr 
2 efi 


k ee a, 

kfjfjCH, 


( 2 ) 


where k ee and k /(/J are electron and muon selection efficiency factors and a accounts for the different trigger 
efficiencies for same-flavour and different-flavour dilepton combinations. The selection efficiency factors 
are calculated using the ratio of dielectron and dimuon events in VRZ according to: 


k 

ee 


kfiM 


A e data (VRZ) 

V data (VRZ)’ 


V data (VRZ) 

V data (VRZ)’ 



(3) 


where e^? g , e^ g and e^ g are the efficiencies of the dielectron, dimuon and electron-muon trigger config¬ 
urations, respectively, and A d ^* /() (VRZ) is the number of ee (pp) data events in VRZ. These selection 
efficiency factors are calculated separately for the cases where both leptons fall within the barrel, both 
fall within the endcap regions, and for barrel-endcap combinations. This is motivated by the fact that the 
trigger efficiencies differ in the central and more forward regions of the detector. This estimate is found to 
be consistent with that resulting from the use of single global k factors, which provides a simpler but less 
precise estimate. In each case the k factors are close to 1.0, and the or estimates obtained using 
k factors from each configuration arc consistent with one another to within 0.2cr. 

The flavour-symmetric background estimate was chosen as the nominal method prior to examining the data 
yields in the signal region, since it relies less heavily on simulation and provides the most precise estimate. 
This data-driven method is cross-checked using the Z boson mass sidebands (mu [81,101] GeV) to fit 
the ??MC events to data in a top control region, CRT. The results are then extrapolated to the signal region 
in the Z boson mass window, as illustrated in Fig. 3. All other backgrounds estimated using the flavour- 
symmetry method are taken directly from MC simulation for this cross-check. Here, Z/y* -t-jets MC events 
are used to model the small residual Z/y* + jets background in the control region, while the jet smearing 
method provides the estimate in the signal region. The normalisation of the tt sample obtained from the fit 
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Figure 3: Diagram indicating the position in the E™ lss versus dilepton invariant mass plane of SR-Z, the control 
region CRT, and the two validation regions (VRT and VRTZ) used to validate the sideband fit for the on-Z search. 
VRT and VRTZ have lower Hj thresholds than CRT and SR-Z. 


is 0.52 ± 0.12 times the nominal MC normalisation, where the uncertainty includes all experimental and 
theoretical sources of uncertainty as discussed in Sect. 7. This result is compatible with observations from 
other ATLAS analyses, which indicate that MC simulation tends to overestimate data in regions dominated 
by tt events accompanied by much jet activity [79,80]. MC simulation has also been seen to overestimate 
contributions from processes in regions with high £™ ss [81]. In selections with high £™ ss but including 
lower Hj, such as those used in the off-Z analysis, this downwards scaling is less dramatic. The results of 
the cross-check using the Z boson mass sidebands are shown in Table 5, with the sideband fit yielding a 
prediction slightly higher than, but consistent with, the flavour-symmetry estimate. This test is repeated 
varying the MC simulation sample used to model the tt background. The nominal Powheg+Pythia tt MC 
sample is replaced with a sample using Alpgen, and the fit is performed again. The same test is performed 
using a Powheg tt MC sample that uses Herwig, rather than Pythia, for the parton shower. In all cases the 
estimates are found to be consistent within lcr. This cross-check using tt MC events is further validated 
in identical regions with intermediate £™ ss (150 < E' | mss < 225 GeV) and slightly looser Hj requirements 
(H t > 500 GeV), as illustrated in Fig. 3. Here the extrapolation in rri(( between the sideband region (VRT) 
and the on-Z region (VRTZ) shows consistent results within approximately lcr between data and the fitted 
prediction. 

The flavour-symmetry method is also tested in these VRs. An overview of the nominal background pre¬ 
dictions, using the flavour-symmetry method, in CRT and these VRs is shown in Fig. 4. This summary 
includes CRT, VRT, VRTZ and two variations of VRT and VRTZ. The first variation, denoted VRT/VRTZ 
(high Hj), shows VRT/ VRTZ with an increased Hj threshold (H\ > 600 GeV), which provides a sample 
of events very close to the SR. The second variation, denoted VRT/VRTZ (high C | mss ), shows VRT/ VRTZ 
with the same £' | mss cut as SR-Z, but the requirement 400 < Hj < 600 GeV is added to provide a sample of 
events very close to the SR. In all cases the data are consistent with the prediction. GGM signal processes 
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Signal 

Flavour-symmetry 

Sideband fit 

region 



SR-Z ee 

2.8 ± 1.4 

4.9+ 1.5 

SR-Z pp 

3.3 ± 1.6 

5.3 + 1.9 


Table 5: The number of events for the flavour-symmetric background estimate in the on-Z signal region (SR-Z) 
using the data-driven method based on data in CRep. This is compared with the prediction for the sum of the 
flavour-symmetric backgrounds (WW, tW, tt and Z —» tt) from a sideband fit to data in CRT. In each case the 
combined statistical and systematic uncertainties are indicated. 


near the boundary of the expected excluded region are expected to contribute little to the normalisation 
regions, with contamination at the level of up to 4 % in CRT and 3 % in VRT. The corresponding con¬ 
tamination in VRTZ is expected to be ~ 10 % across most of the relevant parameter space, increasing to a 
maximum value of ~50 % in the region near m(g) = 700 GeV, p = 200 GeV. 


6.2.2. Flavour-symmetric background in the off-Z search 

The background estimation method of Eq. (2) is extended to allow a prediction of the background dilepton 
mass shape, which is used explicitly to discriminate signal from background in the off-Z search. In 
addition to the k and a correction factors, a third correction factor S (;) is introduced (where i indicates the 
dilepton mass bin): 


Kf{i) = \N^ con (i)k ee aS ee (i), 

N e *(i) = \N^™(i)kwaS^). (4) 

These shape correction factors account for different reconstructed dilepton mass shapes in the ee, pp, and 
ep channels, which result from two effects. First, the offline selection efficiencies for electrons and muons 
depend differently on the lepton p-\ and g. For electrons, the offline selection efficiency increases slowly 
with pt, while it has very little dependence on p \ for muons. Second, the combinations of single-lepton 
and dilepton triggers used for the ee, pp, and ep channels have different efficiencies with respect to the 
offline selection. In particular, for ep events the trigger efficiency with respect to the offline selection at low 
mci is 80%, which is 10-15% lower than the trigger efficiencies in the ee and pp channels. To correct for 
these two effects, tt MC simulation is used. The dilepton mass shape in the ee or pp channel is compared to 
that in the ep channel, after scaling the latter by the a- and k-factor trigger and lepton selection efficiency 
corrections. The ratio of the dilepton mass distributions, N ee {m(()/N efl {mf () or /V /(/ ,( mu)/N e)l (m/■{), is fitted 
with a second-order polynomial, which is then applied as a correction factor, along with a and k, to the 
ep distribution in data. These correction factors have an impact on the predicted background yields of 
approximately a few percent in the ee channel and up to ~ 10— 15 % in the pp channel, depending on the 
signal region. 
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Figure 4: The observed and expected yields in CRT and the VRs in the Z boson mass sidebands (left) and the Z boson 
mass window (right) regions. The bottom plot shows the difference in standard deviations between the observed and 
expected yields. The backgrounds due to WZ, ZZ or rare top processes, as well as from lepton fakes, are included 
under “Other Backgrounds”. 
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The background estimation methodology is validated in a region with exactly two jets and 100 < £ , ' | 1 1ISS < 
150 GeV(VR-offZ). The flavour-symmetric category contributes more than 95 % of the total background 
in this region. The dominant systematic uncertainty on the background prediction is the 6 % uncertainty on 
the trigger efficiency a-factor. The observed dilepton mass shapes arc compared to the SM expectations in 
Fig. 5, indicating consistency between the data and the expected background yields. The observed yields 
and expected backgrounds in the below-Z and above-Z regions are presented in A. For signal models near 
the edge of the sensitivity of this analysis, the contamination from signal events in VR-offZ is less than 
3 %. 



Figure 5: The observed and expected dilepton mass distributions in the electron (left) and muon (right) channel 
of the validation region (VR-offZ) of the off-Z search. Data (black points) are compared to the sum of expected 
backgrounds (solid histograms). The vertical dashed lines indicate the 80 < niff <110 GeV region, which is used 
to normalise the Z + jets background. Example signal models (dashed lines) are overlaid, with m(q), m(x^')/ 
m({)/m(v), and mix j) of each benchmark point being indicated in the figure legend. The bottom plots show the ratio 
of the data to expected background. The error bars indicate the statistical uncertainty in data, while the shaded band 
indicates the total background uncertainty. The last bin contains the overflow. 


6.3. Fake-lepton contribution 

Events from W —> I’v+jets, semileptonic tt and single top ( s- and /-channel) contribute to the background 
in the dilepton channels due to “fake” leptons. These include leptons from b-hadron decays, misidentified 
hadrons or converted photons, and are estimated from data using a matrix method, which is described 
in detail in Ref. [82], This method involves creating a control sample using baseline leptons, thereby 
loosening the lepton isolation and identification requirements and increasing the probability of selecting 
a fake lepton. For each control or signal region, the relevant requirements are applied to this control 
sample, and the number of events with leptons that pass or fail the subsequent signal-lepton requirements 
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are counted. Denoting the number of events passing signal lepton requirements by /V pass and the number 
failing by tVf a n, the number of events containing a fake lepton for a single-lepton selection is given by 

_ N fail - n /e real - 1 )/Vpass 
tV f ak e (i/ e fake _ l/ e real) ’ 

where e fake is the efficiency with which fake leptons passing the baseline lepton selection also pass sig¬ 
nal lepton requirements and e real is the relative identification efficiency (from baseline to signal lepton 
selection) for real leptons. This principle is expanded to a dilepton sample using a four-by-four matrix to 
account for the various possible real-fake combinations for the two leading leptons in the event. 

The efficiency for fake leptons is estimated in control regions enriched with multi-jet events. Events are 
selected if they contain at least one baseline lepton, one signal jet with pj > 60 GeV and low E™ ss 
(<30 GeV). The background due to processes containing prompt leptons, estimated from MC samples, 
is subtracted from the total data contribution in this region. From the resulting data sample the fraction 
of events in which the baseline leptons pass signal lepton requirements gives the fake efficiency. This 
calculation is performed separately for events with ^-tagged jets and those without to take into account the 
various sources from which fake leptons originate. The real-lepton efficiency is estimated using Z —» t + t~ 
events in a data sample enriched with leptonically decaying Z bosons. Both the real-lepton and fake-lepton 
efficiencies are further binned as a function of pj and //. 


6.4. Estimation of other backgrounds 

The remaining background processes, including diboson events with a Z boson decaying to leptons and the 
tt + W( W)/Z and t + Z backgrounds, are estimated from MC simulation. In these cases the most accurate 
theoretical cross sections available are used, as summarised in Table 1. Care is taken to ensure that the 
flavour-symmetric component of these backgrounds (for events where the two leptons do not originate 
from the same Z decay) is not double-counted. 


7. Systematic uncertainties 

Systematic uncertainties have an impact on the predicted signal region yields from the dominant back¬ 
grounds, the fake-lepton estimation, and the yields from backgrounds predicted using simulation alone. 
The expected signal yields are also affected by systematic uncertainties. All sources of systematic uncer¬ 
tainty considered are discussed in the following subsections. 


7.1. Experimental uncertainties 

The experimental uncertainties arise from the modelling of both the signal processes and backgrounds 
estimated using MC simulation. Uncertainties associated with the jet energy scale (JES) are assessed 
using both simulation and in-situ measurements [70,71], The JES uncertainty is influenced by the event 
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topology, flavour composition, jet pj and ?/, as well as by the pile-up. The jet energy resolution (JER) is 
also affected by pile-up, and is estimated using in-situ measurements [83]. An uncertainty associated with 
the JVF requirement for selected jets is also applied by varying the JVF threshold up (0.28) and down 
(0.21) with respect to the nominal value of 0.25. This range of variation is chosen based on a comparison 
of the efficiency of a JVF requirement in dijet events in data and MC simulation. 

To distinguish between heavy-flavour-enriched and light-flavour-enriched event samples, A-jet tagging is 
used. The uncertainties associated with the ^-tagging efficiency and the light/charm quark mis-tag rates 
are measured in tf-enriched samples [84,85] and dijet samples [86], respectively. 

Small uncertainties on the lepton energy scales and momentum resolutions are measured in Z —» ( + {~, 
J/if/ —» l + f~ and W —> ( :± v event samples [64]. These are propagated to the £ , ' | 11ISS uncertainty, along with 
the uncertainties due to the JES and JER. An additional uncertainty on the energy scale of topological 
clusters in the calorimeters not associated with reconstructed objects (the £™ 1SS soft term) is also applied 
to the E™ ss calculation. 

The trigger efficiency is assigned a 5 % uncertainty following studies comparing the efficiency in simula¬ 
tion to that measured in Z —» ( + t~ events in data. 

The data-driven background estimates are subject to uncertainties associated with the methods employed 
and the limited number of events used in their estimation. The Z/y* + jets background estimate has an 
uncertainty to account for differences between pseudo-data and MC events, the choice of seed region 
definition, the statistical precision of the seed region, and the jet response functions used to create the 
pseudo-data. Uncertainties in the flavour-symmetric background estimate include those related to the 
electron and muon selection efficiency factors k ee and k m , the trigger efficiency factor a, and, for the off-Z 
search only, the dilepton mass shape S (i) reweighting factors. Uncertainties attributed to the subtraction 
of the non-flavour-symmetric backgrounds, and those due to limited statistical precision in the ep control 
regions, are also included. Finally, an uncertainty derived from the difference in real-lepton efficiency 
observed in tt and Z —> t + l~ events is assigned to the fake-background prediction. An additional uncer¬ 
tainty due to the number of events in the control samples used to derive the real efficiencies and fake rates 
is assigned to this background, as well as a 20 % uncertainty on the MC background subtraction in the 
control samples. 


7.2. Theoretical uncertainties on background processes 

For all backgrounds estimated from MC simulation, the following theoretical uncertainties are considered. 
The uncertainties due to the choice of factorisation and renormalisation scales are calculated by varying 
the nominal values by a factor of two. Uncertainties on the PDFs are evaluated following the prescription 
recommended by PDF4LHC [87]. Total cross-section uncertainties of 22 % [37] and 50 % are applied to 
tt +W/Z and tt + WW sub-processes, respectively. For the tt + W and tt +Z sub-processes, an additional 
uncertainty is evaluated by comparing samples generated with different numbers of partons, to account 
for the impact of the finite number of partons generated in the nominal samples. For the WZ and ZZ 
diboson samples, a parton shower uncertainty is estimated by comparing samples showered with Pythia 
and Herwig+Jimmy [88,89] and cross-section uncertainties of 5 % and 7 % are applied, respectively. These 
cross-section uncertainties are estimated from variations of the value of the strong coupling constant, the 
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PDF and the generator scales. For the small contribution from t + Z, a 50 % uncertainty is assigned. 
Finally, a statistical uncertainty derived from the finite size of the MC samples used in the background 
estimation process is included. 


7.3. Dominant uncertainties on the background estimates 

The dominant uncertainties in each signal region, along with their values relative to the total background 
expectation, are summarised in Table 6. In all signal regions the largest uncertainty is that associated with 
the flavour-symmetric background. The statistical uncertainty on the flavour-symmetric background due to 
the finite data yields in the e/i CRs is 24 % in the on-Z SR. This statistical uncertainty is also the dominant 
uncertainty for all SRs of the off-Z analysis except for SR-loose, for which the systematic uncertainty 
on the flavour-symmetric background prediction dominates. In SR-Z the combined MC generator and 
parton shower modelling uncertainty on the WZ background (7 %), as well as the uncertainty due to the 
fake-lepton background (14 %), are also important. 


Source 



Relative systematic uncertainty [%] 



SR-Z 

SR-loose 

SR-2j-bveto 

SR-2j-btag 

SR-4j-bveto 

SR-4j-btag 

Total systematic uncertainty 

29 

7.1 

13 

9.3 

30 

15 

Flavour-symmetry statistical 

24 

1.7 

9.3 

6.2 

23 

12 

Flavour-symmetry systematic 

4 

5.7 

6.7 

5.9 

11 

6.6 

Z/y* +jets 

- 

2.1 

6.3 

3.5 

14 

7.0 

Fake lepton 

14 

3.2 

1.4 

1.2 

1.8 

2.2 

WZ MC + parton shower 

7 

- 

- 

- 


- 


Table 6: Overview of the dominant sources of systematic uncertainty on the background estimate in the signal re¬ 
gions. Their relative values with respect to the total background expectation are shown (in %). For the off-Z region, 
the full dilepton mass range is used, and in all cases the ee + fi/j contributions are considered together. 


7.4. Theoretical uncertainties on signal processes 

Signal cross sections are calculated to next-to-leading order in the strong coupling constant, adding the 
resummation of soft gluon emission at NLO+NLL accuracy [55-59]. The nominal cross section and the 
uncertainty are taken from an envelope of cross-section predictions using different PDF sets and factor¬ 
isation and renormalisation scales, as described in Ref. [90]. For the simplified models the uncertainty 
on the initial-state radiation modelling is important in the case of small mass differences during the cas¬ 
cade decays. MadGraph+Pythia samples are used to assess this uncertainty, with the factorisation and 
normalisation scale, the MadGraph parameter used for jet matching, the MadGraph parameter used to set 
the QCD radiation scale and the Pythia parameter responsible for the value of the QCD scale for final- 
state radiation, each being varied up and down by a factor of two. The resulting uncertainty on the signal 
acceptance is up to ~ 25 % in regions with small mass differences within the decay chains. 
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8. Results 


For the on-Z search, the resulting background estimates in the signal regions, along with the observed 
event yields, are displayed in Table 7. The dominant backgrounds are those due to flavour-symmetric 
and WZ and ZZ diboson processes. In the electron and muon channel combined, 10.6 ±3.2 events are 
expected and 29 are observed. For each of these regions, a local probability for the background estimate 
to produce a fluctuation greater than or equal to the excess observed in the data is calculated using pseudo¬ 
experiments. When expressed in terms of the number of standard deviations, this value is referred to as 
the local significance, or simply the significance. These significances are quantified in the last column 
of Table 11 and correspond to a 1.7cr deviation in the muon channel and a 3.0cr deviation in the electron 
channel, with the combined significance, calculated from the sum of the background predictions and 
observed yields in the muon and electron channels, being 3.0cr. The uncertainties on the background 
predictions in the ee and pp channels are correlated as they are dominated by the statistical uncertainty 
of the e/j data sample that is used to derive the flavour-symmetric background in both channels. Since 
this sample is common to both channels, the relative statistical error on the flavour-symmetric background 
estimation does not decrease when combining the ee and pp channels. No excess was reported in the 
CMS analysis of the Z + jets + £ , ' | 1 . 1ISS final state based on yfs = 8 TcV data [24] ; however, the kinematic 
requirements used in that search differ from those used in this paper. 


Channel 

SR-Z ee 

SR-Z pp 

SR-Z same-flavour 



combined 

Observed events 

16 

13 

29 

Expected background events 

4.2+ 1.6 

6.4 ± 2.2 

10.6 + 3.2 

Flavour-symmetric backgrounds 

2.8+ 1.4 

3.3 + 1.6 

6.0 ±2.6 

Z/y* +jets (jet-smearing) 

0.05 ± 0.04 

0 02 +0 03 
u,u -0.02 

0.07 ± 0.05 

Rare top 

0.18 + 0.06 

0.17 + 0.06 

0.35 + 0.12 

WZ/ZZ diboson 

1.2 ±0.5 

1.7 ±0.6 

2.9+ 1.0 

Fake leptons 

O 

1 + 
o o 

1 2 +13 

1 * -1.2 

1 3 +1 - 7 

1 -1.3 


Table 7: Results in the on-Z SRs (SR-Z). The flavour symmetric, Z/y* + jets and fake-lepton background compon¬ 
ents are all derived using data-driven estimates described in the text. All other backgrounds are taken from MC 
simulation. The displayed uncertainties include the statistical and systematic uncertainty components combined. 

Dilepton invariant mass and E™ ss distributions in the electron and muon on-Z SR arc shown in Fig. 6, 
with Hj and jet multiplicity being shown in Fig. 7. For the SR selection a requirement is imposed to reject 
events with AoHjct, 2 , /:. l | mss ) < 0.4 to further suppress the background from Z/y* + jets processes with 
mismeasured jets. 

In Fig. 8, the distribution of events in the on-Z SR as a function of Ar/kjctj 2 , £'™ ss ) (before this requirement 
is applied) is shown. In these figures the shapes of the flavour-symmetric and Z/y* + jets backgrounds are 
derived using MC simulation and the normalisation is taken according to the data driven estimate. 

For the off-Z search, the dilepton mass distributions in the five SRs are presented in Figs. 9 and 10, 
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m M [GeV] m M [GeV] 



Ey iss [GeV] Ej iss [GeV] 

Figure 6: The dilepton mass (top) and ss (bottom) distributions for the electron (left) and muon (right) channel 
in the on-Z SRs after having applied the requirement A</>(jetj 2 ,£™ ss ) > 0.4. All uncertainties are included in the 
hatched uncertainty band. Two example GGM (tan/? = 1.5) signal models are overlaid. For the £,’.' | T " SS distributions, 
the last bin contains the overflow. The backgrounds due to WZ, ZZ or rare top processes, as well as from fake 
leptons, are included under “Other Backgrounds”. The negligible contribution from Z+jets is omitted from these 
distributions. 
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Figure 7: The Hj (top) and jet multiplicity (bottom) distributions for the electron (left) and muon (right) channel 
in the on-Z SRs after having applied the requirement A</>(jetj 2 , £.j nlss ) > 0.4. All uncertainties are included in the 
hatched uncertainty band. Two example GGM (tan/? = 1.5) signal models are overlaid. For the H\ distributions, 
the last bin contains the overflow. The backgrounds due to WZ, ZZ or rare top processes, as well as from fake 
leptons, are included under “Other Backgrounds”. The negligible contribution from Z+jets is omitted from these 
distributions. 
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Figure 8: The distribution of the A <p between the leading jet and £j mss (top) and the sub-leading jet and /t’.j r " ss 
(bottom) for the electron (left) and muon (right) channel in the on-Z SRs before having applied the require¬ 
ment A0(jetj 2 ,£™ ss ) > 0.4. All uncertainties are included in the hatched uncertainty band. Two example GGM 
(tan /j = 1.5) signal models are overlaid. The backgrounds due to WZ, ZZ or rare top processes, as well as from fake 
leptons, are included under “Other Backgrounds”. 
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and summarised in Fig. 11. The expected backgrounds and observed yields in the below-Z and above-Z 
regions for SR-2j-bveto, SR-4j-bveto, and SR-loose are presented in Tables 8, 9, and 10, respectively. 
Corresponding results for SR-2j-btag and SR-4j-btag are presented in A. The data are consistent with the 
expected SM backgrounds in all regions. In the SR-loose region with 20 < ma < 70 GeV, similar to the 
region in which the CMS Collaboration observed a 2.6cr excess, 1133 events are observed, compared to 
an expectation of 1190 ± 40 ± 70 events. 
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Figure 9: The observed and expected dilepton mass distributions in the off-Z SR-loose (top), SR-2j-bveto (middle), 
and SR-4j-bveto (bottom). The vertical dashed lines indicate the 80 < <110 GeV region, which is used to 

normalise the Z + jets background and is thus not treated as a search region. Example signal models (dashed lines) 
are overlaid, with m(q)/m(g), /m(j(f), m(t)/m(v), 28d m(x®) of each benchmark point being indicated in the 

figure legend. The last bin contains the overflow. All uncertainties are included in the hatched uncertainty band. 
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Figure 10: The observed and expected dilepton mass distributions in the SR-2j-btag (top) and SR-4j-btag (bottom) 
signal regions of the off-Z search. The vertical dashed lines indicate the 80 < ma <110 GeV region, which is 
used to normalise the Z + jets background and is thus not treated as a search region. Example signal models of 
squark- or gluino-pair production (dashed lines) are overlaid, with m(g), m(jfy/m(x f), m(£)/m(v), and m(x °) of 
each benchmark point being indicated in the figure legend. The last bin contains the overflow. All uncertainties are 
included in the hatched uncertainty band. 
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Figure 11: The observed and expected yields in the below-Z (left) and above-Z (right) dilepton mass regions, for the 
VR and five SRs of the off-Z search. Here below-Z is 20 < mu < 70 GeV for VR-offZ and SR-loose and otherwise 
20 < mu < 80 GeV; above-Z is mu >110 GeV. The bottom plot shows the difference in standard deviations 
between the observed and expected yields. Results are shown for the ee and pp channels as well as for the sum. 


9. Interpretation of results 

In this section, exclusion limits are shown for the SUSY models described in Sect. 3. The asymptotic 
CLs prescription [91], implemented in the HistFitter program [92], is used to determine upper limits 
at 95 % confidence level (CL). All signal and background uncertainties are taken into account using 
a Gaussian model of nuisance parameter integration. All uncertainties except that on the signal cross 
section arc included in the limit-setting configuration. The impact of varying the signal cross sections by 
their uncertainties is indicated separately. Numbers quoted in the text are evaluated from the observed 
exclusion limit based on the nominal signal cross section minus its lcr theoretical uncertainty. 

For the on-Z analysis, the data exceeds the background expectations in the ee (pp) channel with a signi¬ 
ficance of 3.0 (1.7) standard deviations. Exclusion limits in specific models allow us to illustrate which 
regions of the model parameter space are affected by the observed excess, by comparing the expected and 
observed limits. The results in SR-Z ee and SR-Z pp (Table 7) are considered simultaneously. The signal 
contamination in CRe/u is found to be at the ~ 1 % level, and is therefore neglected in this procedure. The 
expected and observed exclusion contours, in the plane of p versus m{g) for the GGM model, are shown in 
Fig. 12. The ±lo- exp and ±2cr eX p experimental uncertainty bands indicate the impact on the expected limit 
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Below-Z (20 < m {e < 80 GeV) 

SR-2j-bveto ee 

SR-2j-bveto /j/j 

SR-2j-bveto same-flavour 
combined 

Observed events 

30 

24 

54 

Expected background events 

26 + 4 + 3 

24 + 4 + 3 

50 + 8 + 5 

Flavour-symmetric backgrounds 

Z/y* + jets 

Rare top 

VV’Z/ZZ diboson 

Fake leptons 

24 + 4 + 3 

0.6 + 0.3 + 0.7 
< 0.1 

0.6+ 0.2+ 0.1 

0.6+ 0.9+ 0.1 

22 + 4 + 3 

1.6+ 0.6 + 1.4 
<0.1 

0.6 + 0.2 + 0.2 
<0.1 

46 + 8 + 4 
2.2+ 0.7+ 1.7 
< 0.1 
1.2+ 0.3+ 0.2 
0.2 + 0.9 + 0.1 


Above-Z (mu >110 GeV) 

SR-2j-bveto ee 

SR-2j-bveto /j/j 

SR-2j-bveto same-flavour 
combined 

Observed events 

26 

29 

55 

Expected background events 

35 + 5 ± 4 

38+4 + 8 

73 + 9 + 9 

Flavour-symmetric backgrounds 

Z/y* + jets 

Rare top 

VV'Z/ZZ diboson 

Fake leptons 

33+4 + 4 

0.3 ± 0.2 + 0.3 
< 0.1 

0.3 + 0.1 +0.1 
1.7+ 1.1 +0.2 

30 + 4 + 3 

5.6+ 0.6+ 7.5 
<0.1 

0.6+ 0.2+ 0.1 

1.3+ 1.1 +0.5 

63 + 8 + 5 
5.9 + 0.7 + 7.5 
< 0.1 

0.8+ 0.2+ 0.1 
3.0 ± 1.5 + 0.4 


Table 8: Results in the off-Z search region SR-2j-bveto, in the below-Z range (20 < mu < 80 GeV, top) and above-Z 
range (mu >110 GeV, bottom). The flavour symmetric, Z/y* + jets and fake lepton background components are 
all derived using data-driven estimates described in the text. All other backgrounds are taken from MC simulation. 
The first uncertainty is statistical and the second is systematic. 


of all uncertainties considered on the background processes. The ±lcr t S h ^ r y uncertainty lines around the 
observed limit illustrate the change in the observed limit as the nominal signal cross section is scaled up 
and down by the theoretical cross-section uncertainty. Given the observed excess of events with respect to 
the SM prediction, the observed limits are weaker than expected. In the case of the tan /)’ =1.5 exclusion 
contour, the on-Z analysis is able to exclude gluino masses up to 850 GeV for fi > 450 GeV, whereas 
gluino masses of up to 820 GeV arc excluded for the tan fi = 30 model for p > 600 GeV. The lower 
exclusion reach for the tan [i — 30 models is due to the fact that the branching fraction for jp® —> ZG is 
significantly smaller at tan /j = 30 than at tan fi =1.5. 

For the off-Z search, the limits for the squark-pair (gluino-pair) model are based on the results of SR-2j- 
bveto (SR-4j-bveto). The yields in the combined ee+/j/j channels are used. Signal contamination in the 
e/j control region used for the flavour-symmetry method is taken into account by subtracting the expected 
increase in the background prediction from the signal yields. For each point in the signal model parameter 
space, limits on the signal strength arc calculated using a “sliding window” approach. The binning in 
SR-2j-bveto (SR-4j-bveto) defines 45 (21) possible dilepton mass windows to use for the squark-pair 
(gluino-pair) model interpretation, of which the ten (nine) windows with the best expected sensitivity arc 
selected. For each point in the signal model parameter space, the dilepton mass window with the best 
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Below-Z (20 < m u < 80 GeV) 

SR-4j-bveto ee 

SR-4j-bveto pp 

SR-4j-bveto same-flavour 
combined 

Observed events 

1 

5 

6 

Expected background events 

4.7+ 1.6 ± 1.1 

3.6+ 1.5 + 1.0 

8.2 ±3.1 ± 1.4 

[•favour-symmetric backgrounds 
Z/y* + jets 

Rare top 

VV'Z/ZZ diboson 

Fake leptons 

4.1 ± 1.6+ 1.1 
0.4 ± 0.2 ± 0.3 
< 0.1 
< 0.1 

0.2 ± 0.3 ± 0.0 

3.5 ± 1.5 ± 1.0 

0.0 ± 0.0 ± 0.4 
<0.1 
<0.1 
<0.1 

7.7 ±3.1 ± 1.3 
0.4 ± 0.2 ± 0.5 
< 0.1 
< 0.1 

0.1 ±0.3 ±0.0 


Above-Z (mu >110 GeV) 

SR-4j-bveto ee 

SR-4j-bveto pp 

SR-4j-bveto same-flavour 
combined 

Observed events 

2 

9 

11 

Expected background events 

5.7 ± 1.6+ 1.2 

4.5 ±1.3 ±1.7 

10 + 3 + 2 

Flavour-symmetric backgrounds 
Z/y* + jets 

Rare top 

WZ/ZZ diboson 

Fake leptons 

5.5+ 1.6+ 1.2 
0.2 ±0.1 ±0.1 
< 0.1 
< 0.1 
< 0.2 

4.3+ 1.3 + 1.0 

0.0 ±0.0+ 1.3 
<0.1 

0.2 ±0.1 ±0.0 
<0.1 

9.8 ±2.9+ 1.4 
0.2 ±0.1 ± 1.3 
< 0.1 

0.2 ±0.1 ±0.0 
< 0.2 


Table 9: Results in the off-Z search region SR-4j-bveto, in the below-Z range (20 < mu < 80 GeV, top) and above-Z 
range (mu >110 GeV, bottom). Details are the same as in Table 8. 


expected limit on the signal strength is selected. The excluded regions in the squark-LSP and gluino-LSP 
mass planes are shown in Fig. 13. The analysis probes squarks with masses up to 780 GeV, and gluinos 
with masses up to 1170 GeV. 

The signal regions in these analyses are also used to place upper limits on the allowed number of BSM 
events (/Vbsm) in each region. The observed (5^ s ) and expected (S^p) 95 % CL upper limits are also 
derived using the CL$ procedure. These upper limits on /Vbsm can be interpreted as upper limits on the 
visible BSM cross section ((ecr)^? ) by normalising Vbsm by the total integrated luminosity. Here (ecr)^ s 
is defined as the product of the signal production cross section, acceptance and reconstruction efficiency. 
The results are obtained using asymptotic formulae [93] in the case of the off-Z numbers. For SR-Z, 
with a considerably smaller sample size, pseudo-experiments are used. These numbers are presented in 
Table 11 for the on-Z search. Model-independent upper limits on the visible BSM cross section in the 
below-Z and above-Z ranges of the five signal regions in the off-Z search are presented in Tables 12 and 
13, respectively. Limits for the most sensitive dilepton mass windows of SR-2j-bveto and SR-4j-bveto 
used for the squark- and gluino-pair model interpretations, respectively, are presented in Tables 14 and 15. 
These tables also present the confidence level observed for the background-only hypothesis CLg, and the 
one-sided discovery p-value, p(s = 0), which is the probability that the event yield obtained in a single 
hypothetical background-only experiment (signal, s = 0) is greater than that observed in this dataset. The 
p(s = 0) value is truncated at 0.5. 
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Below-Z (20 < ma < 70 GeV) 

SR-loose ee 

SR-loose pp 

SR-loose same-flavour 
combined 

Observed events 

509 

624 

1133 

Expected background events 

510 + 20 + 40 

680 ±20 + 50 

1190 + 40 + 70 

Flavour-symmetric backgrounds 

Z/y* + jets 

Rare top 

WZ/ZZ 

Fake leptons 

490 ±20 + 40 

2.5 ±0.8 ±3.2 

0.3 ± 0.0 ± 0.0 

1.1 ±0.3 ±0.1 

16 + 4 + 2 

650 ±20 + 50 

8 + 2 + 5 

0.4 ± 0.0 ± 0.0 

1.2 ±0.2 ±0.4 

23 + 5 + 1 

1140 ±40 + 70 

11+2 + 7 

0.7 ± 0.0 ± 0.0 
2.4 ± 0.4 ± 0.4 

38 + 6 + 4 


Above-Z (niff >110 GeV) 

SR-loose ee 

SR-loose pp 

SR-loose same-flavour 
combined 

Observed events 

746 

859 

1605 

Expected background events 

760 ±20 + 60 

830 ±20 + 70 

1600 ± 40 ± 100 

Flavour-symmetric backgrounds 

Z/y* + jets 

Rare top 

WZ/ZZ diboson 

Fake leptons 

730 ±20 + 60 

0.9 ±0.2+ 1.1 

0.2 ± 0.0 ± 0.0 

0.6 ± 0.2 ± 0.2 

30 + 5 + 5 

800 ±20 + 60 
21+3 + 24 

0.2 ± 0.0 ± 0.0 

1.0 ±0.2 ±0.1 

6.7 ± 3.7 ± 1.7 

1500 + 40+ 100 
22 + 3 + 24 

0.4 ±0.0 ±0.0 
1.6 ±0.3 ±0.2 
37 + 6 + 5 


Table 10: Results in the off-Z search region SR-loose, in the below-Z range (20 < niff < 70 GeV, top) and above-Z 
range (niff >110 GeV, bottom). Details are the same as in Table 8. 


Signal region 

Channel 

<«r>& M 

S 95 h 

obs 

c95 
^ exp 

CL b 

II 

o 

Gaussian significance 

SR-Z 

ee + pp 

1.46 

29.6 

12 

0.998 

0.0013 

3.0 


ee 

1.00 

20.2 

8+ 2 

0.998 

0.0013 

3.0 


pp 

0.72 

14.7 

9-2 

0.951 

0.0430 

1.7 


Table 11: From left to right: 95 % CL upper limits on the visible cross section ((ecr) 9 ? ) and on the number of 
signal events (S 9 jj s ); the expected 95 % CL upper limit on the number of signal events is denoted by S jjj, and is 
derived from the expected number of background events (and the ±lcr uncertainty on the expectation); two-sided 
CLb value, which is the confidence level observed for the background-only hypothesis; the discovery p-value for 0 
signal strength s (p(s = 0)), and the Gaussian significance for the on-Z search. 
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GGM: higgsino-like tan|5 = 1.5, M = M 2 = 1 TeV, m(q)=1.5 TeV 




Figure 12: The 95 % CL exclusion limit from the on-Z combined same-flavour channels in the /r versus m(g) plane 
in the GGM model with tan/1 =1.5 (top) and tan/? - 30 (bottom). The dark blue dashed line indicates the expected 
limits at 95 % CL and the green (yellow) bands show the + lcr (+2<x) variation on the expected limit as a consequence 
of the experimental and theoretical uncertainties on the background prediction. The observed limits are shown by 
the solid red lines, with the dotted red lines indicating the limit obtained upon varying the signal cross section by 
+ 1 <x. The region below the grey line has the gluino mass less than the lightest neutralino mass and is hence not 
considered. The value of the lightest neutralino mass is indicated by the x-axis inset. 
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Figure 13: Excluded region in the (top) squark-LSP mass plane using the SR-2j-bveto results and (bottom) gluino- 
LSP mass plane using the SR-4j-bveto results. The observed, expected, and ±lcr expected exclusion contours are 
indicated. The observed limits obtained upon varying the signal cross section by ±lcr are also indicated. The region 
to the left of the diagonal dashed line has the squark mass less than the LSP mass and is hence not considered. Three 
signal benchmark points are shown, with their SUSY particle masses indicated in parentheses. 
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Signal region 

Channel 

^Vdata 

^V”bkg 


c 95 
°obs 

c 95 
^ exp 

CL b 

o 

II 

SR-2j-bveto 

ee + pp 

54 

50 ± 8 ± 5 

1.38 

28.0 

24^ 

0.66 

0.35 


ee 

30 

26 + 4 + 3 

0.99 

20.1 

18 + ^ 

0.73 

0.28 


pp 

24 

24 + 3 + 3 

0.88 

17.8 

18-6 

0.50 

0.50 

SR-2j-btag 

ee + pp 

79 

104+ 11+7 

0.98 

19.8 

30 + g 0 

0.06 

0.50 


ee 

40 

49 + 6 + 4 

0.85 

17.2 

20i 

0.19 

0.50 


W 

39 

56 + 6 + 5 

0.63 

12.8 

20^ 

0.06 

0.50 

SR-4j-bveto 

ee + pp 

6 

8.2 ±3.1 ± 1.4 

0.38 

7.7 


0.37 

0.50 


ee 

1 

4.7 + 1.6+ 1.1 

0.19 

3.9 

5-4!?;? 

0.08 

0.50 


pp 

5 

3.6+ 1.5 ± 1.0 

0.41 

8.4 

6 5 +2 - 9 

u.j_l l 

0.75 

0.26 

SR-4j-btag 

ee + pp 

31 

38 + 6 + 3 

0.85 

17.3 

19 + _l 

0.25 

0.50 


ee 

14 

18 + 3 + 2 

0.51 

10.3 

13+5 

15 i 

0.30 

0.50 


pp 

17 

20 + 4 + 2 

0.54 

10.9 

0.33 

0.50 

SR-loose 

ee + pp 

1133 

1190 + 40 + 70 

6.82 

138.4 

170-40 

100^0 

HO-30 

0.28 

0.50 


ee 

509 

510 + 20 + 40 

4.88 

99.0 

0.51 

0.48 


W 

624 

680 ±20 + 50 

4.10 

83.3 

0.18 

0.50 


Table 12: Summary of model-independent upper limits for the five signal regions, in the below-Z region (20 < 
mu < 70 GeV for SR-loose, 20 < mu < 80 GeV for all other signal regions), in the combined ee + pp and individual 
ee and pp channels. Left to right: the observed yield (AWa), total expected background (AW), 95 % CL upper 
limits on the visible cross section ((ecr)^ s ) and on the number of signal events (S^ ). The fifth column (,S'^ p ) 
shows the 95 % CL upper limit on the number of signal events, given the expected number (and ±lcr excursions 
on the expectation) of background events. The last two columns indicate the CLg value, i.e. the confidence level 
observed for the background-only hypothesis, and the discovery p-value (pis - 0)). For an observed number of 
events lower than expected, the discovery p-value has been truncated at 0.5. 


36 



Signal region 

Channel 

^data 

-^bkg 


obs 

c 95 
^ exp 

CL b 

o 

ll 

SR-2j-bveto 

ee + /i/i 

55 

73 ± 9 ± 9 

0.96 

19.4 

27 + y 

0.11 

0.50 


ee 

26 

35 ± 5 ± 4 

0.60 

12.1 

isi 

0.14 

0.50 


W 

29 

38 ±4 ±8 

0.89 

18.1 

20!« 

0.24 

0.50 

SR-2j-btag 

ee + ii[i 

164 

171 ±14+16 

2.19 

44.4 

48^2 

0.39 

0.50 


ee 

83 

81 + 7 + 7 

1.45 

29.5 

28.3^n 0 

0.56 

0.43 


W 

81 

90 + 7 + 14 

1.49 

30.2 

36^g° 

0.33 

0.50 

SR-4j-bveto 

ee + [ i[i 

11 

10 + 3 + 2 

0.56 

11.4 

10 -3 

0.61 

0.42 


ee 

2 

5.7 + 1.6+ 1.2 

0.20 

4.1 

6-0-1.8 

0.13 

0.50 


m 

9 

4.5 ± 1.3 ± 1.7 

0.61 

12.3 

7.7^2 

0.91 

0.08 

SR-4j-btag 

ee + [i[i 

41 

36 + 6 + 5 

1.27 

25.7 

20! 9 

0.72 

0.29 


ee 

23 

18 + 3 + 2 

0.96 

19.5 

15^ 

0.83 

0.17 


m 

18 

19 + 3 + 4 

0.85 

17.2 

17 3 

0.50 

0.50 

SR-loose 

ee + [i/i 

1605 

1600 ± 40 ± 100 

10.58 

214.8 

210 +30 

zlu -40 

0.62 

0.40 


ee 

746 

760 ±20 + 60 

6.63 

134.6 

140+5° 

lw -40 

0.42 

0.50 


[i[i 

859 

830 ±20 + 70 

8.23 

167.1 

150 +so 

0.64 

0.32 


Table 13: Summary of model-independent upper limits for the five signal regions, in the above-Z (w« >110 GeV) 
dilepton mass range, in the combined ee+/i/i and individual ee and /./// channels. Details are the same as in Table 12. 


m u range [GeV] 

^data 

^bkg 

C^obJ^l 

c 95 
^ obs 

a 95 

J exp 

cl b 

pis = 0) 

<' 385 

jy 665,265 

yy745,25 

Slg 

20-50 

35 

26 ± 6 ± 3 

1.32 

26.9 

20!’ 

0.85 

0.15 

17.1 ± 1.6 

3.7 ± 0.4 

0.6 ±0.1 

20-80 

54 

50 ± 8 ± 4 

1.38 

28.0 

24*1 

0.66 

0.35 

[38.0 ±2.4] 

10.4 ± 0.6 

2.1 ±0.2 

50-80 

19 

23 ± 5 ± 2 

0.63 

12.8 

17*| 

0.30 

0.50 

20.9 ± 1.8 

6.7 ±0.5 

1.5 ±0.2 

50-140 

34 

46 ± 7 ± 6 

0.83 

16.9 

20!* 

0.14 

0.50 

27.3 ± 2.0 

28.5 ± 1.0 

6.9 ±0.3 

50-200 

51 

75 ± 9 ± 8 

0.89 

18.1 

26!* 

0.05 

0.50 

28.2 ±2.1 

50.6 ± 1.4 

14.2 ± 0.5 

110-200 

32 

52 ± 7 ± 7 

0.69 

14.1 

20!* 

0.05 

0.50 

2.8 ±0.6 

[34.0 ± 1.1] 

10.5 ± 0.4 

170-260 

12 

24 ± 5 ± 2 

0.40 

8.2 

12*1 

0.03 

0.50 

0.4 ±0.2 

14.8 ± 0.7 

11.9 ±0.4 

170-290 

16 

26 ± 5 ± 2 

0.43 

8.7 

13*_l 

0.08 

0.50 

0.4 ±0.2 

16.1 ±0.8 

16.8 ±0.5 

>170 

25 

34 ± 6 ± 3 

0.68 

13.9 

19!* 

0.15 

0.50 

0.4 ±0.2 

18.5 ±0.8 

[25.7 ± 0.6] 

>230 

16 

13.1 ±3.2 ±2.3 

0.88 

17.9 

14*1 

0.72 

0.29 

0.3 ±0.2 

5.0 ±0.4 

17.8 ±0.5 


Table 14: Summary of model-independent upper limits for SR-2j-bveto, in the combined ee + [ip and individual ee 
and pp channels, for the 10 dilepton mass windows used for the squark-pair interpretation. Details are the same as 
in Table 12. The last three columns indicate the expected signal yield for three squark-pair model benchmark points; 
the first (second) number indicates the squark (LSP) mass. The signal yield in square brackets indicates the best 
selected dilepton mass window for the given benchmark point. 
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^data 

Nbk% 


c95 
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o95 
^ exp 

cl b 

O 

II 

<>> 

yy825,585 

sig 

^1025,545 

iV sig 

pr 1185,65 
sig 

20-50 

4 

3.1 ±2.3 ±0.9 

0.40 

8.2 

7.5 + _f° 

0.70 

0.38 

4.4 ±0.7 

0.8 ± 0.1 

0.1 ± 0.0 

20-80 

6 

8.2 ± 3.1 ± 1.4 

0.38 

7.7 

8.3T« 

0.37 

0.50 

[ 12.8 ± 1 . 1 ] 

2.0 ± 0.2 

0.2 ± 0.0 

50-140 

6 

8.2 ± 2.7 ± 1.4 

0.37 

7.5 

8 . 2 ^ 

0.35 

0.50 

21.4 ± 1.4 

4.9 ±0.3 

0.7 ±0.1 

110-200 

9 

5.6 ± 2.3 ± 1.4 

0.59 

12.0 

o 4 + 3.5 
-2.0 

0.85 

0.17 

4.2 ±0.6 

6.3 ± 0.3 

1.0 ± 0.1 

140-260 

6 

5.0 ±2.1 ±0.8 

0.43 

8.6 

7 4 + 30 
'■^-1.4 

0.66 

0.38 

1.3 ±0.4 

[8.0 ±0.4] 

1.6 ± 0.1 

>20 

17 

18 ± 4 ± 3 

0.63 

12.8 

+ 1 

xf- 

0.46 

0.50 

27.4 ± 1.6 

14.4 ±0.5 

7.4 ±0.2 

>140 

7 

7.2 ± 2.4 ± 1.3 

0.41 

8.3 

8 - 2 ^i i 

0.52 

0.50 

1.6 ±0.4 

8.6 ±0.4 

6.7 ± 0.2 

>200 

2 

4.8 ± 1.8 ± 1.1 

0.21 

4.2 

5.9!“ 

0.23 

0.50 

0.4 ±0.2 

4.2 ±0.3 

6.0 ± 0.2 

>260 

1 

2.3 ±1.2 ±0.7 

0.19 

3.9 

4 2 +L9 

^• z -0.3 

0.34 

0.50 

0.3 ± 0.2 

0.7 ±0.1 

[5.1 ±0.1] 


Table 15: Summary of model-independent upper limits for SR-4j-bveto, in the combined ee + /j/j and individual ee 
and /.(/j channels, for the nine dilepton mass windows used for the gluino-pair interpretation. Details are the same 
as in Table 12. The last three columns indicate the expected signal yield for three gluino-pair model benchmark 
points; the first (second) number indicates the gluino (LSP) mass. The signal yield in square brackets indicates the 
best selected dilepton mass window for the given benchmark point. 
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10. Summary 


This paper presents results of two searches for supersymmetric particles in events with two same-flavour 
opposite-sign leptons, jets, and £™ ss , using 20.3 fb' 1 of 8 TeV pp collisions recorded by the ATLAS 
detector at the LHC. The first search targets events with a lepton pair with invariant mass consistent 
with that of the Z boson and hence probes models in which the lepton pair is produced from the decay 
Z —> It. In this search 6.4 ± 2.2 (4.2 ± 1.6) events from SM processes are expected in the pp ( ee ) SR-Z, as 
predicted using almost exclusively data-driven methods. The background estimates for the major and most 
difficult-to-model backgrounds are cross-checked using MC simulation normalised in data control regions, 
providing further confidence in the SR prediction. Following this assessment of the expected background 
contribution to the SR the number of events in data is higher than anticipated, with 13 observed in SR-Z pp 
and 16 in SR-Z ee. This corresponding significances are 1.7 standard deviations in the muon channel and 
3.0 standard deviations in the electron channel. These results are interpreted in a supersymmetric model of 
general gauge mediation, and probe gluino masses up to 900 GeV. The second search targets events with a 
lepton pair with invariant mass inconsistent with Z boson decay, and probes models with the decay chain 
X® I n this case the data are found to be consistent with the expected SM backgrounds. No 

evidence for an excess is observed in the region similar to that in which CMS reported a 2.6cr excess [24]. 
The results are interpreted in simplified models with squark- and gluino-pair production, and probe squark 
(gluino) masses up to about 780 (1170) GeV. 
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A. Additional results of off-Z search 


This section provides additional results of the off-Z search. The expected backgrounds and observed yields 
in the below-Z and above-Z regions for VR, SR-2j-btag, and SR-4j-btag, are presented in Tables 16, 17, 
and 18, respectively. 


Below-Z (20 < m„ < 70 GeV) 

VR-offZ ee 

VR-offZ /i/i 

VR-offZ same-flavour 
combined 

Observed events 

465 

742 

1207 

Expected background events 

445 ± 15 ± 36 

682 ±23 + 53 

1128 + 37 + 69 

Flavour-symmetric backgrounds 

Z/y* + jets 

Rare top 

VV'Z/ZZ diboson 

Fake leptons 

425 + 15 + 36 

1.5 ±0.5 + 1.6 

0.1 ±0.0 ±0.0 

0.6 ±0.2 ±0.1 
18+4 + 2 

661 +22 + 53 

2.7 ±0.6 ±3.5 

0.1 ±0.0 ±0.0 

0.8 ±0.2 ±0.1 

18 + 4 + 4 

1086 ±37 + 68 
4.1 ±0.8 ±4.7 

0.1 ±0.0 ±0.0 
1.4 ±0.2 ±0.2 

36 ± 5 ± 7 


Above-Z (mu >110 GeV) 

VR-offZ ee 

VR-offZ mi 

VR-offZ same-flavour 
combined 

Observed events 

550 

732 

1282 

Expected background events 

594+ 18+48 

696 + 21 +55 

1290 ±38+79 

Flavour-symmetric backgrounds 

Z/y* + jets 

Rare top 

WZ/ZZ diboson 

Fake leptons 

571 ±17 + 48 

1.9 ±0.7 ±2.0 
< 0.1 

0.9 ± 0.2 ± 0.2 

21+4 + 2 

684 + 21 +55 

3.8 ±0.4 ±6.0 

< 0.1 

0.9 ± 0.2 ± 0.2 

7.9 ±3.1 ±2.9 

1254 + 38 + 79 

5.7 ±0.8 ±7.5 

< 0.1 

1.8 ±0.3 ±0.2 

29 ± 5 ± 4 


Table 16: Results in the off-Z validation region (VR-offZ), in the below-Z range (20 < mu < 70 GeV, top) and 
above-Z range (mu >110 GeV, bottom). The flavour symmetric, Z/y* + jets and fake lepton background com¬ 
ponents are all derived using data-driven estimates described in the text. All other backgrounds are taken from MC 
simulation. The first uncertainty is statistical and the second is systematic. 
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Below-Z (20 < mu < 80 GeV) 

SR-2j-btag ee 

SR-2j-btag /./// 

SR-2j-btag same-flavour 
combined 

Observed events 

40 

39 

79 

Expected background events 

49 + 6 + 4 

56 ± 6 ± 5 

104+ 11+7 

Flavour-symmetric backgrounds 

Z/r* + jets 

Rare top 

W'Z/ZZ diboson 

Fake leptons 

45 + 5+4 

1.8 ± 1.0 ±0.8 

0.1 ±0.0 ±0.0 
< 0.1 
2.3+ 1.2 ±0.3 

49 + 6 + 5 

3.1 ± 1.3+ 1.9 

0.1 ±0.0 ±0.0 

0.1 ±0.1 ±0.1 

3.4+ 1.9 ±0.2 

94+11+7 
4.9+ 1.6 ±2.2 

0.1 ±0.0 ±0.0 
0.1 ±0.1 ±0.1 
5.7 ±2.3 ±0.6 


Above-Z {mu >110 GeV) 

SR-2j-btag ee 

SR-2j-btag /./// 

SR-2j-btag same-flavour 
combined 

Observed events 

83 

81 

164 

Expected background events 

81 + 7 + 7 

90 + 7+14 

171 ±14+16 

Flavour-symmetric backgrounds 

Z/y* + jets 

Rare top 

VV’Z/ZZ diboson 

Fake leptons 

78 + 7 + 7 

0.8 ± 0.5 ± 0.4 
< 0.1 
< 0.1 
2.4+ 1.6 ±0.8 

77 + 7 + 7 

11 + 1 + 13 
< 0.1 
< 0.1 

1.6+ 1.3 ±0.2 

155+13 + 10 
12+1 + 13 

0.1 ±0.0 ±0.0 
< 0.1 
4.0 ±2.1 ±0.7 


Table 17: Results in the off-Z search region SR-2j-btag, in the below-Z range (20 < < 80 GeV, top) and above-Z 

range {mu >110 GeV, bottom). Details are the same as in Table 8. 
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Below-Z (20 < mu < 80 GeV) 

SR-4j-btag ee 

SR-4j-btag yU/u 

SR-4j-btag same-flavour 
combined 

Observed events 

14 

17 

31 

Expected background events 

18 + 3 + 2 

20 + 4 + 2 

38 + 6 + 3 

Flavour-symmetric backgrounds 

Z/r* + jets 

Rare top 

W'Z/ZZ diboson 

Fake leptons 

17 + 3 + 2 

0.5 ± 0.3 ± 0.5 
< 0.1 
< 0.1 

0.3 ± 0.6 ± 0.0 

18 + 3 + 2 

0.7 ± 0.3 ± 0.9 
< 0.1 
< 0.1 

1.3+ 1.2 ±0.0 

35 + 6 + 3 
1.2 ±0.4+ 1.1 

0.1 ±0.0 ±0.0 
< 0.1 
1.6+ 1.4 ±0.2 


Above-Z (mu >110 GeV) 

SR-4j-btag ee 

SR-4j-btag /j/.i 

SR-4j-btag same-flavour 
combined 

Observed events 

23 

18 

41 

Expected background events 

18 + 3 + 2 

19 + 3 + 4 

36 + 6 + 5 

Flavour-symmetric backgrounds 

Z/y* + jets 

Rare top 

VV’Z/ZZ diboson 

Fake leptons 

17 + 3 + 2 

0.2 ±0.1 ±0.2 
< 0.1 
< 0.1 

< 0.6 

16 + 3 + 2 

2.4 ± 0.3 ± 4.0 
< 0.1 
< 0.1 

0.3 ±0.6 ±0.1 

33 + 6 + 3 
2.7 ±0.3 ±4.1 
< 0.1 
< 0.1 

0.0 ± 0.9 ± 0.2 


Table 18: Results in the off-Z search region SR-4j-tag, in the below-Z range (20 < nin < 80 GeV, top) and above-Z 
range (mu >110 GeV, bottom). Details are the same as in Table 8. 
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